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ABSTRACT 



Aims. We search for evidence of a relation between properties of young stellar objects (YSOs) and their parent molecular clouds to 
understand the initial conditions of high-mass star formation. 

Methods. A sample of 135 sources was selected from the Infrared Astronomical Satellite (IRAS) Point Source Catalog, on the ba- 
sis of their red color to enhance the possibility of discovering young sources. Using the Kolner Observatorium fiir SubMillimeter 
Astronomie (KOSMA) 3-m telescope, a single-point survey in '^CO 7 = 2-1 was carried out for the entire sample, and 14 sources 
were mapped further. Archival mid-infrared (MIR) data were compared with the '""CO emissions to identify evolutionary stages of 
the sources. A "CO observed sample was assembled to investigate the correlation between '^CO line width of the clouds and the 
luminosity of the associated YSOs. 

Results. We identified 98 sources suitable for star formation analyses for which relevant parameters were calculated. We detected 18 
cores from 14 mapped sources, which were identified with eight pre-UC Hii regions and one UC Hn region, two high-mass cores 
earlier than pre-UC Hii phase, four possible star forming clusters, and three sourceless cores. By compiling a large (360 sources) 
'""CO observed sample, a good correlation was found between the "CO line width of the clouds and the bolometric luminosity of the 
associated YSOs, which can be fitted as a power law, Ig (AVb/ km s"') = (-0.023 ± 0.044) + (0.135 ± 0.012) Ig (Lboi/Lo). Results 
show that luminous (> 10' Lq) YSOs tend to be associated with both more massive and more turbulent (AV^ > 2 km s"') molecular 
cloud structures. 

Key words, stars: formation — ISM: clouds — ISM: molecules — ISM: kinematic and dynamics 



1. Introduction 

The past decade has witnessed significant progress in the study 
of high-mass star formatio n. Observations at millimeter and 
submillimeter wavelengths ('Zhang et al.', 1998L iBeuther et al.' 



2002; Keto, 2002; Zhang, 2005; Cesaroni^etaLl |2007) suggest 
that massive proto B stars can form by disk mediated accre- 
• tion, which is similar to the scenario that produces low-mass 
' stars. However, most of the studies focus on relatively evolved 
stages, when the central star has already formed and hydrogen 
burning has begun, characterized by surrounding ultra compact 

i UC) H II regions a nd strong emission from complex molecules 
Churchwell |2002|) . In contrast, the extremely early stages are 
poorly understood to date. In particular, knowledge to evolution- 
ary stages prior to the onset of H n regions are crucial to under- 
standing the initial conditions of high-mass star formation. 

It is known that stars are formed in molecular clouds. 
Therefore, the relation between forming stars and parent clouds 
is important to understand the formation process and the proper- 
ties of the eventual stars. On galaxy scales, star formation activ- 
ities are usually described by the so-called Schmidt law, which 
relates the star formation rate (SFR) to the su rface density of 
gas: SsFR Sf^,^, where the index = 1 - 2 dSchmidd Il959t: 
^ennicutt, 1998; Gao & Solomon, 200^. Studies of Galactic 
dense cores have shown that this relation ma y be universal an d 
can be connected to Galactic star formation dWu et all l2005ah . 



ILarso ^ (Il981h studied the turbulence in star forming clouds and 
found a strong correlation between the internal velocity disper- 
sion cr of the region and its size L: cr(kms"') ocL(pc)'*^**. This 
relation, also called the Larson law, is valid for low-mass cores 
but is found to be break down in high-mass core s > lO^M p 
(Casein & Myers, 1995; Plume etal., 1997; Guan etalt |2008|). 
This is indicative of the different status of turbulence in low- and 
high-mass cores. The breakdown of the Larson law can be inter- 
preted as evidence of widespread supersonic turbulence in high- 
mas s cores, i n con trast to subsonic turbulent low-mass cores 
(Plu me et al.Lll99'7 r). A molecular line width is an observational 
indicator of turbulence in clouds, and bolometric luminosity is 
an indicator of forming stars. Any relation between these quan- 
tities may help us to understand the initial star forming process. 

Here we report results from a '''CO J -2-1 survey towards 
135 IRAS sources using the KOSMA 3-m telescope. To search 
for high-mass star forming regions in their early stages, we se- 
lect a sample on the basis of their red IRAS color to enhance the 
possibility of finding young sources. We present the primary re- 
sults and investigate the relation between line width in molecular 
clouds and bolometric luminosity of associated infrared sources. 
We describe our sample selection in Sect. |2] and observations in 
Sect.|3] In Sect.|4]we present statistical results of the single-point 
survey (Sect. |4.1] i and follow-up mapping (Sect. |4.2] i. We discuss 
the Ay - L relation as well as other relations in Sect.|5] and sum- 
marize the paper in Sect.|6] 
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2. Sample 

We selected the sample from the Inf rared Astronomical Sat ellite 
(IRAS) Point Source Catalog (PSC, ' Beichman et"aD 'l988') ver- 
sion 2.1 according to our developed color criteria (Wuetal], 
l2003b . namely: 

(a) /lOO/j«, < 500 Jy, \g{f25^m I flliim) > 0.7, 

lg(/60^,m//i2/im) > 1-4, whcrc is the flux density; 

(b) lack of 6 cm radio continuum radiation to exclude 
potential Hii associations; 

(c) declination 5 > -20°, so that targets are accessible to the 
telescope KOSMA. 

Criterion (a) was chosen so that the sample sources were red- 
der and possibly fain ter, thus may be younger tha n those selected 
based on traditional IWood & Churchwelll d 1989k) color criteria. 
Criterion (b) helps to exclude any known Hn regions brighter 
than current detection limit. Therefore, the sample should rep- 
resent extremely young stellar objects (YSOs), mostly at evo- 
lutionary stages earlier than the UC Hii phase. The 6 cm 
radio continuum data was extracted from three surveys: (1) 
0° < ^ < 75° 4.85 GHz r adio continuum survey completed 
by [Gregory & Condoril (119911) with the 91-m NRAO telescope; 
(2) -29° < S < 9.5 ° 4.85 GHz radio continuum survey led by 
iGriffith et all (11994 with the 64-m Parkes telescope; and (3) 
-9.5° < ^ < 10° 4.85 GHz radio continuum survey led by 
iGriffith et all (1 19951) wifli the 64-m Parkes telescope. 

Criteria (a) and (c) lead to 500 sources being selected from 
the PSC, which contains 245,889 sources. However, only 135 
sources were observed because of limited observing time and 
after applying criterion (b). These sources represent the sam- 
ple reported in this paper. The sample sources are concentrated 
across the Galactic plane and cover a wide range of longitude, 
10° < / < 230°. 



3. Observations 

A single-point survey in '^^CO 7 = 2-1 (220.398 GHz) was car- 
ried out from September 2002 to March 2003 using the Kolner 
Observatorium fiir SubMillimeter Astronomic (KOSMAQ) 3-m 
telescope on Gornergrat near Zermatt in Switzerland. All of the 
sample sources were surveyed in '^CO 7 = 2-1. About half 
of the sample sources were also observed in '^CO 7 = 2-1 
(230.538GHz) and 14 of them were mapped in '^CO 7 = 2-1. 

The beamwidth of the KOSMA at 230 GHz was 130". 
The pointing accuracy was superior to 10". The telescope was 
equipped with a dual-channel SIS receiver, which had a noise 
temperature of 150 K. A high resolution spectrometer with 
2048 channels was employed and the spectral resolution was 
165.5 KHz, giving a velocity resolution of 0.22 km s"^ The main 
beam temperature (r^b) had been corrected for the effects of 
Earth's atmosphere, antenna cover loss, radiation loss, and for- 
ward spillover and scattering efficiency (92%). From the cali- 
brated Jupiter observations, the main beam efficiency ri„,i, was 
estimated as 68% during our observation. On-the-fly mode was 
adopted during mapping, with a mapping step of 60". Most 
maps were extended until the line intensity decreased to half 
of the maximum value or even lower The GILDAS0 software 



' The KOSMA 3 m radiotelescope at Gornergrat-Siid Observatory 
is operated by the University of Cologne and supported by special 
funding from the Land NRW. The Observatory is administered by the 
Internationale Stiftung Hochalpine Forschungsstationen Jungfraujoch 
und Gornergrat, Bern. 

^ available at http : //www . iram . f r/IRAMFR/GILDAS | 
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Fig. 1: Example spectra of '-^CO 7 = 2-1 towards the IRAS 
sources given in the text. 



package (CLASS/GREG/SIC) was used for the data reduction 
(Guilloteau & Lucas, .2000) . 



4. Results 

4.1. Survey 

Among the entire sample of 135 IRAS sources, we identified 98 
sources suitable for star formation analyses (another 37 sources 
were excluded either because they had multiple components 
or bad baselines, or failed to be detected), of which 60 have 
both '^CO 7 = 2-1 and '^CO 7 = 2-1 data. Figure [1] presents 
example spectra of '^CO 7 = 2-1: (a) IRAS 00117+6412, a 
perfect Gaussian profile; (b) IRAS 02541+6208, a fairly nar- 
row line width; (c) IRAS 06067+21 38, a broad line with red 
wing, also seen in 7 = 1 - transition (Wue t al.Ll2003b ; (d) IRAS 
20326+3757, a blue wing; (e) IRAS 18278-0212, red asymme- 
try; (f) IRAS 21379+5106, two peaks; (g) IRAS 19348+2229, 
two components; and (h) IRAS 02485+6902, multiple compo- 
nents. 

Observed and derived parameters are listed in an online 
Table[T] starting with IRAS name and its J2000 equatorial coordi- 
nates in Cols. (1) to (3). By Gaussian fit, we obtain the observed 
parameters including main beam temperature T^b, local stan- 
dard of rest velocity Vlsris, and '^CO 7 = 2-1 line width (full 
width at half-maximum) AVu for each source, listed in Cols. (4) 
to (7). When a line profile is obviously non-Gaussian, the pa- 
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rameters are measured with a cursor (e.g., IWu & Evansll2003h . 
and the velocity uncertainty is given as the velocity resolution; 
when the line profile has distinctive multiple components, only 
the strongest component is shown, indicated by a character m in 
corresponding Tmb columns. 

The distance to most sources was unavailable in the litera- 
ture. The kinematic distances were calculated based on the radial 
vel ocity Vlsrh and the v elocity field of the outer Galaxy given 
bv lBrand& Bliti (Il993h . When two kinematic distances were 
available, we selected the closer one, except when the closer dis- 
tance is too small (<100pc). For 8 sources, however, no reason- 
able distances could be calculated in this way and we assumed 
that the distance to these sources is 1 kpc. These are marked as * 
in the distance Col. (8) of Table[T] 

The bolometric luminosity was calculated based on the dis- 
tances and the IRAS fluxes in four bands (12, 25, 60, 100 /im), 
following the formula given bv lCasoli et al.l ( ll986h 

Lboi = 5.4d2(/i2^„,/0.79 + f25^J2 + f6o^„/3.9 + fm^,n/9.9) Lo, 

where D is the distance in kpc and is the flux density in Jansky. 
The uncertainties in luminosity originate in the kinematic dis- 
tances and the quality of the IRAS source fluxes. Most of the 
sample sources have high or moderate quality in all the four 
bands. Twenty-one sources with upper limit fluxes in one or two 
bands are marked as u luminosity in Col. (9) of Table[T] 

Assuming local thermodynamic equilibrium (LTE) and that 
the '^CO 7 = 2-1 transition is optically thin (i.e. T13 < 1), we 
derive excitation temperatures, optical depth and column den- 
sities for '^CO, using radiation transfer equation dGarden et al.L 
1199 li) . Based on the assumption of LTE, '^CO and '^CO share 
the same excitation temperature Tex, which can be derived from 
the main beam temperature of optically thick '^CO, Tmbia- When 
Ti3 > 1, an optical depth correction factor Cj - Tn/{1 - e"'"") 
is multiplied by its corresponding column density. The rel- 
ative CO abundance ['^CO/Ha] is estim ated to extend from 
2.5 X 10-5 (iRodriguez et al.LfT982h to lO'^^ (iGarden et al.L[T99Tl) . 
and we adopt the median value of 6.25 x 10"^. Using the ter- 
restrial [i^C/i^C] ratio of 89, we adopt a value for ['^CO/Hj] of 
7.0x10"^ when computing the column density of H2. These pa- 
rameters are listed in Cols. (10) to (13). References of former 
works are given in the last Col. (14) of Table[T] 

The distribution of '^CO 7 = 2-1 line width of this sam- 
ple has a mean of 3.09 km s ' and a standard deviation of 
1.06kms This line width is relatively smaller than that of 
typical brig ht/red IRAS sou rces associated with water masers 
(3.5kms~^ IWu et alj|200ll : note that this value was measured 
in 7 = 1 - transition), while significantly larger than that of a 
molecular clou d hosting interm ediate-mass star formation activi- 
ties (~2kms~ MSun et al.l2006l averaged throughout the Perseus 
cloud). The luminosities are distributed over a wide range, from 
20 Lq to about 10^ Lq, with a mean of 10"* L©. The high disper- 
sion of luminosities indicates that these sources are embedded 
in very different environments. This luminosity distributio n is 
similar to the young 'low' sources of iMolinari et aljri996l (see 
their Fig. 6), in agreement with the assumption that our sam- 
ple group may be relatively younger than that chosen by tradi- 
tional color criteria. The excitation temperature Tgx ranges from 
4.4 to 22.5 K, with an average of 9.7 K. This suggests that very 
cold gases surround the sample source s, colder than thos e sur- 
rounding the luminous IRAS sources (IZhu & WuL I200I . The 
'^CO column densities are (1.2 - 28.7) xlO'^cm"^, with an av- 
erage of 6.2 xlO'5 cm"^, while H2 column densities are (1.7 - 
40.8) xlO^' cm"^, with an average of 8.9 xlO^' cm"^. These den- 



sities are roughly close to th e critical value for gravitational col- 
lapse dHartquist et al.Lll998h . 

4.2. Mapping 

To improve our understanding of the properties of the sur- 
veyed sample, 14 sources were mapped in '^CO 7 = 2-1 and 
compared with archival mid-infrared (MIR) continuum data. 
Mapped sources were selected from the surveyed sample as 
those with only single emission component, and they almost 
evenly cover longitude 70° < I < 230°, avoiding low Galactic 
longitudes, where '^CO lines are often affected by multiple ve- 
locity components from the Galactic molecular ring. Using these 
sources as a guide, maps were extended until at least one core 
was resolved. We name a map on the basis of its guide source 
name, as outlined in Fig.|2] In four cases, one map resolved two 
cores, resulting in 18 cores in total. We found that 13 cores are 
associated with the original guide sources, two cores are associ- 
ated with other IRAS sources, and three cores have no embedded 
infrared source (sourceless hereinafter). A core is named after its 
associated IRAS source; for a sourceless core, it is named after 
its nearest IRAS source plus relative direction to the core (e.g., 
20067+3415NE). See Table|2]for core properties. 

The core size (Col. 2 of Table |2]i is defined as an equivalent 
linear size R = y/A/n, where A is the projected area of each 
cloud within the 50% contour (highlighted in Fig. |2|. It is cor- 
rected for the effect of beam smearing by multiplying its value 

by a factor -^^^bs ~ ^mb/^^s' where flobs is the angular diameter 
of the core and flmb is the beamwidth. For three cores, the ob- 
served angular diameters are comparable to the beamwidth, so 
that the cores are just marginally resolved and the correspond- 
ing core sizes are highly uncertain. In a few cases, maps were 
not complete to 50% of the peak intensity, and can only infer 
lower limits to R (indicated by a symbol '>'). The average line 
width of each core (Col. 3) is determined by combining all the 
spectra in the core and then fitting a Gaussian profile to the av- 
erage spectrum. In a few cases, the average spectra show line 
asymmetry/absorption and need to be fitted with two Gaussian 
profiles, and then the line width of the stronger profile is given. 
The typical uncertainty in the average line width is 0.04 km s"'. 
Column (4) lists the luminosity also given in Table [T] for ref- 
erence. Peak volume densities for H2, n(H2) (Col. 5), and the 
LTE core masses, Mlte (Col. 6), are calculated based on both R 
and the peak '^CO column densities determined by interpolat- 
ing the maps. For three maps (IRAS 06067 -h2138, 07024-1 102, 
and 21391-H5802), however, no A^('^CO) are available in Table[l] 
because of a lack of '^CO data. To estimate their core prop- 
erties, we assume reasonable excitation temperatures: for IRAS 
06067 and 07024, we assume a typical Tgx of 15 K; and for IRAS 
21391, we assume that Tex equals the dust temperature (25 K, 
iBeltran et al.ll2002l) . Column (7) present s the virial mass derive d 
from the sizes and line widths following lMacLaren et al.l ( fl988l) . 
The ratio of virial to LTE mass a - MyiJMuE is listed in the last 
Col. (8) of Table |2] We exclude the marginally resolved cores 
when computing averages except for the line width column. 

Overall, the core mass ranges from ~ 10-Mo to IO'^Mq, the 
linear size from 0. 1 1 pc to 2.41 pc, and molecular hydrogen den- 
sity is in the range ~ 10^ - lO"* cm"^. The luminosities are once 
again, distributed across a wide range, from 30 Lq to 1.7xlO"^L0. 
Overall, the line width AVi^ > ~ 2 km s and has an average of 
2.80kms"', smaller than that of the entire surveyed sample. We 
find an average value of 1.3 for the ratio of virial to LTE core 
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Table 2: Core Properties 



Core 


R 


AVi3 


Lbo\ 


Peak n(H2) 


Mlte 




a' 




(pc) 


(kms"') 


(lO'Lo) 


(10' cm-3) 


(Me) 


(Mo) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


00117+6412 


0.57 


2.58 


4.29 


5.13 


1.8E+2 


4.8E+2 


2.7 


00557+5612 


1.08 


2.05 


1.47" 


1.66 


6.6E+2 


5.7E+2 


0.9 


03101+5821 


>1.08 


2.52 


1.20 


1.76 


>5.8E+2 


8.6E+2 


1.5 


03260+31 IP 


>0.23 


3.05 


0.29 


15.96 


>1.6E+2 


2.6E+2 


1.6 


03260+31 UNE^ 


0.11 


2.51 




15.96 


8.0E+1 


8.8E+1 


1.1 


03414+3200 


>0.34 


1.92 


0.05" 


5.35 


>8.5E+1 


1.6E+2 


1.8 


05168+3634 


>2.41 


2.86 


17.13 


1.57 


>1.2E+4 


2.4E+3 


0.2 


05168+3634SW^ 


0.60" 


2.61 




1.57 


3.0E+3 


5.2E+2 


0.2 


06067+2138 


0.55 


3.36'' 


0.03" 


3.83"^ 


2.5E+2 


7.8E+2 


3.1 


06103+1523 


1.97 


2.56 


9.49 


1.02 


2.8E+3 


1.6E+3 


0.6 


07024-1102 


>o.ir" 


1.99 


0.57 


64.73"^ 


>9.0E+0 


5.8E+1 


6.4 


20067+3415 


1.76 


3.80'' 


1.14 


1.14 


3.9E+3 


3.2E+3 


0.8 


20067+3415N& 


0.45™ 


3.67 




1.14 


9.7E+2 


7.5E+2 


0.8 


20149+3913 


0.43 


3.19 


0.30" 


2.62 


1.4E+3 


5.5E+2 


0.4 


20151+3911" 


>1.30 


3.83 


0.70" 


2.62 


>4.1E+3 


2.4E+3 


0.6 


21391+5802 


0.32 


2.78 


0.26 


11.87"^ 


1.3E+2 


3.1E+2 


2.4 


22198+6336 


1.22 


2.2^ 


1.53" 


2.39 


1.7E+3 


7.5E+2 


0.4 


22506+5944 


1.31 


2.88 


6.83 


1.40 


l.OE+3 


1.4E+3 


1.4 


Average^ 


0.98 


2.80 


3.02 


4.95 


1.9E+3 


l.lE+3 


1.3 



Note. — " not guide IR AS source. Lumino sity calculated based on IRAS fluxes and distance: 03260+3111 at 0.4 kpc jHariu et^kl |199^ and 
20151+3911 at 1.7kpc jMotte et alll2007h . '' average spectrum fitted with two Gaussian profiles, list the stronger one. no '^CO data, density 
calculated by assuming reasonable T^x (see text). ™ core's angular diameter comparable to beamwidth, marginally resolved. sourceless core. 
" upper limit. ^ average calculating does not include values of marginally resolved cores, except for the line width column. * a = MyiJMuE- 



mass, a. Overall the mapping sample infers a ~ I, indicating 
that most of the cores appear to be virialized. 

Comparisons between '^CO maps and MIR images are pre- 
sented in Fig.|2l and a detailed evolutionary identification of in- 
dividual mapped sources is presented in the online Appendix A. 

5. Discussion 

5.1. The line width-luminosity relation 

The empirical correlation of fine width versus luminosity has 
been found by va r ious author s from observations in C"*0 
Saito et al.', 2001; Ridge etal., 200 3b, as well as in NH^ 
(Wouterloot et al., 1988rMvers et al.. ll99U iHarju et all 1 19931 



iLadd et al.L ll994; Jijina et al., 1999D. The fine width is in general 
found to increase with lu minosity, for differe nt lg{AV)/ Ig(Lboi) 
slopes: 0.13-0.19 f or NH3 jJiiina et al.L[T999t) and 0. 1 1 for C'^O 
jSaito et al.L 12001). Given the relatively wide availability of the 
archival '■'CO data, it is helpful to compile an up-to-date '^CO 
observed sample to investigate the luminosity-line width relation 
in case of '^CO. 

In Fig. [3] we plot in logarithmic space the line width 
versus luminosity from our sample an d other ' ^CO observed 



samples a dopted fro m the literature dBeichman et alj 



Dent et all 1985; Fis cher et all Il985t lYamashita et al] 



1986; 



1989; 



Wu et alil2001i.i2004tiRidge et al.Ll200l . This sample contains 



360 sources in total. One finds that the luminosity of the IRAS 
sources is well correlated with the '^CO line width, as fitted by 
a power law: 



^^'^ ' (-0.023 ±0.044) + (0.135 ± 0.012) Ig ( — 



I kms 



where the correlation coefficient c.c. - 0.69. This suggests that 
the mass of the forming stellar objects is linked to the dynamic 



status of their parent clouds. [Saito et alj ( 1200 ll) measured a sim- 
ilar correlation in the Centaurus tangential region and suggested 
that the mass of the formed stars is determined by the internal 
velocity dispersion of the dense cores. If the velocity dispersion 
is a reliable indicator of the turbulence, this is consistent with 
the idea that turbulence is different in high-mass and low-mass 
cores. We note that the '^CO line widths adopted from the litera- 
ture were measured from 7=1-0 transition with smaller beams. 

It is generally agreed that embedded infrared point sources 
of high luminosity (> 10^ Lq) but without associated H 11 regions 
are good ca ndidates to be h igh-mass YSOs in the pre-UC Hii 
phase (e.g., I Wu et al.l |2006| and references therein). However, 
sufficiently young massive objects are not necessarily bright 
at infrared wavelengths; some of them have no infrared coun- 
terparts. High-mass objects at sufficiently young evolutionary 
stages could be quite faint in infrared ranges either because they 
are not yet mature enough to have developed infrared emis- 
sion or they are em bedded very deeply in cold dust. For in- 
stance, ISridharan et alj (120051) identified several 1 .2 mm emit- 
ting high-mass starless cores (HMSCs) that exhibit absorption 
or no emission at the MIR wavelengths; a centimeter-emitting 
UC H ii region was also found without an infrared counterpart 
d^brichetal., 2008). In our mapped sample, faint (< IO-'Lq) 
sources associated with very massive (> IO^Mq) cores (20149, 
20 1 5 1 ) do exist. Although we cannot rule out the possibility that 
these two sources could evolve to only low-mass stars, it is very 
likely that the clouds will eventually fragment to form high-mass 
stars, given the large amount of gas therein. 

Hereafter, for clarity, our luminosity criterion refers to bolo- 
metric luminosity Lboi ^ IO-'Lq, and our line-width criterion 
refers to line width Ay('-'CO y = 2-l)> 2 kms"'. According 
to the Ay - L relation above, luminosity Lboi - 10^ cor- 
responds to AVn - 2.42;';g4j kms"'. Because high-mass stars 
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Fig. 2: '■'CO J -2-1 integrated intensity contours overlay the Midcourse Space Experiment (MSX: IFrice et al.ll200lh band A 
(8.28 fim) images as background, if available. Maps are listed in order of name, except the two without MSX data. IRAS 03258+3 104 
has no band A data and we use band C instead; IRAS 00557+5612 and 22198+6336 both have no MSX data available, as labeled 
on the top of their relevant sub-figures. IRAS point sources are denoted by squares while MSX point sources by triangles. Small 
crosses represent observed points, while large crosses denote the original guide sources (Sect. 14. 2]) . as labeled on the top of each 
map. Dash lines schematically separate resolved cores. For integration, only intensity over three times of the standard deviation 
(3a-) is considered. Contour levels begin from 30% to 90% by 10% of the peak intensity, while 50% level is highlighted by a solid 
thick contour. For MSX images, the grey scale wedge is shown on the right side; the unit is W m"^ sr"' . 
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Fig. 3: Line width plotted versus bolometric luminosity of 360 '^CO observed sources: surveyed sources (open circles), mapped 
sources (filled circles), and other samples adopted from the literature (diamonds). Solid line represents a least squares fitting to the 
data, and the dash lines represent the luminosity/line-width criteria (Sect. lS.ll) . 



are far more luminous than their low-mass counterparts, lumi- 
nous IRAS sources (> 10^ L©) are likely to be high-mass stel- 
lar objects. Therefore, we tentatively (not exclusively) suggest 
the lower limit, AVn - 2kms"', as a characteristic value for 
the line-width criterion, analogous to the widely used luminosity 
criterion. Objects with line width larger than this characteristic 
value are probably high-mass objects. The line-width criterion 
includes 94.5% sources that also satisfy the luminosity criterion 
(Fig. [3]), which implies that line width may be a key parameter in 
measuring the masses of the forming stellar objects in the cores, 
at least in our sample. We note that this criterion (2kms"') is 
larger tha n the typical line width of low mass cores (1.3kms 
iMvers et a l. 1983), and sma ller than the av erage line width of 
high-mass cores (3.5kms''. IWuet aLllMl . 

Applying the luminosity/line-width criteria to our sample, 
there are 68 sources satisfying luminosity criterion, 65 (95.6%) 
of which also satisfy line-width criterion. We suggest that the 
65 sources are candidate high-mass star formation regions in a 
pre-UC H ii phase. For the remaining 30 less luminous sources, 
23 of them satisfy the line-width criterion but not the luminosity 
criterion. We suggest that the 23 sources are high-mass YSO 
candidates earlier than pre-UC H ii phase. 

5.2. Core masses and line widths 

The core masses provide a direct test of our line-width criterion. 
Figure |3] includes 15 cores with luminosities listed in Table |2] 
We exclude sourceless cores and marginally resolved cores in 
the discussion in this section, because the former's luminosity 
cannot be determined and latter's estimated size and mass are 
highly uncertain. The remaining 14 luminosity-available cores 
can be divided into two groups: Group I, which do not satisfy 
the luminosity criterion, including 03414, 21391, 03260, 20149, 
06067, 20151; and group II, satisfying the luminosity criterion, 
including 00557, 22198, 03101, 06103, 00117, 05168, 22506, 
20067 (in order of increasing line width). All group II cores 



also satisfy the line-width criterion, and they are located in the 
upper right panel of Fig. [3] They are very massive, their es- 
timated masses being higher than several IO^Mq, except core 
00117 (1.8 X IO^Mq). On the other hand, group I cores mostly 
satisfy the line-width criterion. They are relatively less massive 
than group II cores, their masses being typically ~ lO^M© or 
more, with two very massive cores 20149 and 20151 (> IO^'Mq). 
The only case that does not satisfy the line-width criterion, core 
03414, has the lowest mass in group I. We note that all group 
I core s are still significa ntly more massive than the low-mass 
cores (iMvers et al.lll983h . 

The high-mass nature of group I cores confirms again (in ad- 
dition to the previously mentioned cores 20149, 20151) that the 
luminosity criterion cannot be applied to some young sources. 
On the other hand, the line-width criterion is applicable to our 
sample. While in terms of inferring mass, line width may not be 
as direct an indicator as luminosity, it is helpful when luminos- 
ity is unavailable or is affected by large uncertainty (e.g., due to 
distance ambiguity, flux upper limit), which is often the case. In 
addition, line width can be measured observationally more easily 
and more accurately than luminosity. 

In Fig.|4](a), we plot lg(Ayi3) versus IgMLXE- A weak cor- 
relation is evident in the data, with a correlation coefficient of 
0.38. This may indicate that, for molecular clouds with associ- 
ated YSOs, the '^^CO line width at some degree is related to the 
cloud mass, and massive cores tend to have larger line widths. 
This weak correlation, together with the strong AV - L corre- 
lation, indicates that massive stars are more likely to form in 
massive molecular cores. The core mass and associated IRAS 
luminosity in our mapped sample are indeed well correlated 
(c.c. - 0.76). However, this correlation may need to be corrected 
for distance effects (0.3-6.08 kpc for mapped sample) because 
both mass and luminosity are proportional to D^. Nevertheless, 
strong mass-luminosity correlations were reported in other re- 
gions or sam ples that have far smaller distance differences than 
our sample tPobashi et al.LII996l:[Saito et al.Ll200ltlRidge etsli 
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Fig. 4: Relations between '■'CO line width AVo and other phys- 
ical parameters of molecular cores: (a) LTE mass Mlte, and (b) 
size R. In both cases, the filed circles represent data and solid 
lines are linear fitting results to the data. The fitting function and 
correlation coefficient are labeled in each panel. 



l2003h . Figure. |4] (b) plots line width AV^ and core size R in 
logarithm. With an average virial mass of 1.1 X 10^ Mo, the 
cores exhibit no correlation between size and line width. This 
indicates that the Larson law is invalid for o ur mapped sample, 
consistent with the results of previous works dPlume et al.lll997l 
(Mvir) > 3.8 X 10^ Mq; Guan et al. 2001 <Mvir) = 5.6 x 10^ Mq). 
The correlations in Fig.|4]are unaffected by distance, because the 
line width and distance of the plotted cores are not correlated 
(c.c. = 0.04). 

We conclude that, based on the currently available sample, 
YSOs with higher bolometric luminosity (> lO"' L©), tend to be 
associated with more massive molecular cloud structures, which 
are usually more turbulent, and have a large '''CO line width, 
AVo > 2kms"'. It is important to note that, the characteristic 
value (2kms"') may not be universal, and can vary from re- 
gion to region and/or from line to line. Further mapping of more 
clouds, as well as higher angular resolution data if available, are 
required to examine the line-width criterion proposed here. 

6. Summary 

We have carried out a '^CO 7 = 2-1 survey of 135 IRAS sources 
selected as potential YSOs earlier than UC Hii regions. Our 
main findings are summarized as follows: 

1. Ninety-eight sources have good enough emission pro- 
file for analysis; some of them show asymmetric line profiles 



of CO 7 = 2-1. The line width is 3.09 km s" and excitation 
temperature is 9.7 K, on average. The H2 column densities are 
(1.7 - 40.8) xlO^' cm"^. Sixty-five sources are suggested to be 
candidate precursors of UC H 11 regions. 

2. Fourteen sources were mapped and resolved as eighteen 
cores, which have been identified with eight pre-UC H 11 regions 
and one UC H 11 region, two high-mass cores earlier than pre-UC 
H II phase, four possible star forming clusters, and three source- 
less cores. 

3. For molecular clouds with known associated YSOs and 
measured Lbob '"'CO line width AVo of the clouds is corre- 
lated with the bolometric luminosity of the YSOs. Based on 
the current '^'CO observed sample (360 sources in total), this 



correlation can be fitted as a power law, lg(Ayi3/kms ') = 
(-0.023 + 0.044) + (0.135 + 0.012) lg(Lboi/L0). 

4. Luminous (> 10^ L©) YSOs tend to be produced in more 
massive and more turbulent (AVn > 2kms"') molecular cloud 
structures. 

5. High-mass stars are more likely to form in massive molec- 
ular clouds. 
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Appendix A: Individual Analyses 

In Fig.|2] we compare the integrated intensity of '^CO emissions 
to IRAS (point sources) and MSX (point sources and images) 
data. We use the most sensitive band (band A, centered on 
8.28 jum) images of MSX when available. We note that maps in 
Fig.|2]are labeled by the original guide sources (Sect. |4.2] l. IRAS 
03258+3104 has no band A data and we use band C (centered 
on 12.13yum) instead. IRAS 00557+5612 and 22198+6336 
both have no MSX data. Individual analyses of each map are 
presented as follows. 

IRAS 00117+6412: '^CO emission peak coincides well 
with IRAS and MSX point sources, and there are also strong 
counterparts in all four MSX bands. This distinctive '^CO core 
is massive (1.8 x IQ^Mq), in agreement with the conclusion 
deduced from the luminosity/line-width criteria (Sect. 15. it . We 
therefore suggest that it is a pre-UC Hii region. Strong 22 
GHz water maser dWouterlo ot et al., 1993) and outflow activity 
jZhang et al.Ll2005t Zhao et al..,2003i) have been detected within 
this area, providing evidence of active star formation. 

IRAS 00557+5612: MSX data are unavailable close to 
this region, but the IRAS source matches well to the ^^CO 
core peak. A core as massive as 6.6 x lO^M© agrees with the 
conclusion deduced from the luminosity/line-width criteria. 
We therefore suggest that it is a pre-UC H ii region. A velocity 
gradient of 0.36km s"' pc"' from northeast to southwest is 
inferred, yielding a rotating angular velocity of 1.16 x 10'*'* s~ 



Accor ding to '^CO 7=1-0 and HCO+ mapping bv lZhu & Wul 
(12001) . two subcores exist within this core. 

IRAS 03101+5821: '^CO emission coincides with infrared 
point sources and image as well. A core more massive than 
5.8 X IO^Mq agrees with the conclusion drawn from the 
luminosity/line-width criteria. We therefore suggest that it is a 
pre-UC H ii regi on. A 22 GHz water m aser has been detected 
within this area dWouterloot etaUll993h . 

IRAS 03258+3104: a MSX band A image is unavailable for 
this region and we use band C instead. '^CO emission within 
this area is more diffuse than that in former sources. At least 
two cores (03260+3111 and 03260+31 1 INE) are resolved 
within an area of 0.34pc. The larger core coincides with IRAS 
03260+31 1 1, which does not satisfy our color selection criteria 
and w as suggested as an UC Hii region by Churchwell et alj 
(Il990l) . Taking both their line width and luminosity into account, 
we suggest that core 03260+3111 is a high-mass object in UC 
Hii phase, while 03260+311 INE is a sourceless core. The 
original guide source of the map, IRAS 03258+3104, is not 
associated with any resolved ^^CO core. It has been suggested 
to be a Class object driving a low-mass bipolar CO outflow 
(iKnee & SandellLl2000l) . 

IRAS 03414+3200: ^^CO is quite diffuse across the entire 
area of 0.34pc, so that no distinctive '■'CO core is found. 
However, several infrared point sources are evident close to 
the 90% contour within one beam, superimposed on a steep 
density gradient. Although its mass (> 85 Mq) and line width 
1.92kms"') are relatively low in all the mapped sources, it 
appears to be a star forming cluster 

IRAS 05168+3634: at least two cores (05168+3634 and 
05168+3634SW) are present within 3.01 pc. The dominant 



northeastern core appears to be associated with the infrared 
sources. The total core mass higher than 1.5x IO^'Mq agrees with 
the conclusion deduced from the luminosity/line-width criteria. 
We suggest that the dominant core is a high-mass star forming 
region in pre-UC Hii phase, and the SW core is a sourceless 
core. This is consistent with results from Molinari et al. ( 199 
Strong outflow activity was identified (Br and et ali 1 199' 
Zhang et al., 2005) at the position of the NE core, and the 
outflow driving source appears deviated to the infrared source 
IRAS 05168+363 4. A 22 GHz water maser was detected by 
iPalla etalJ(ll991h in this region. 

IRAS 06067+2138: '^CO core coincides with the IRAS 
point source but is without MSX counterpart. The core mass 
Mlte (2.5 X IO^Mq) is only one third of its virial mass Myir 
(7.8 X IO^Mq), indicating that the core is not yet gravitationally 
bound. Taking its large line width (3.36 km s"') into account, 
we suggest that it is a high-mass object earlier than pre-UC H ii 
phase. 

IRAS 06103+1523: '^CO emission coincides with both 
infrared point sources and image. A core as massive as 
2.8 X 10-^ Mo agrees with the conclusion deduced from the 
luminosity /line-width criteria. IRAS 06103+1523 is found to be 
two point sources by MSX data, implying that a fine structure 
may exist there. A denser molecular tracer (e.g., N2H^ or 
HCO^) and a higher resolution (several arcsec) are needed to 
study these fine structures. We therefore suggest that it is a 
high-mass star forming cluster. 

IRAS 07024-1102: '^CO map is incomplete but a core is 
clearly evident. The core coincides with the IRAS point source 
but does not have a MSX counterpart. Although its luminosity 
(570 Lq) is a little lower than the luminosity criterion, it does 
generate a line width (1.99 km s"') very close to the line-width 
criterion. We therefore suggest that it is a high-mass object 
earlier than pre-UC H n phase. 

IRAS 20067+3415: the '^CO gas distribution is quite 
complex in this region. At least two cores (20067+3415 and 
20067 +34 15NE) are revealed. The dominant southwestern core 
coincides with infrared point sources, while the northeastern 
core is a sourceless core. Several sub-structures are revealed 
within an area of 2.21 pc with total mass of 4.9 x IO^Mq. Thus, 
we suggest that core 20067+3415 is a high-mass star forming 
cluster The MIR emission and luminosity are relatively weak 
compared to those of other pre-UC Hii regions, indicating a 
very early evolutionary stage. 

IRAS 20149+3913: '^CO emission reveals two cores 
(20149+3913 and 20151+3911); both also coincide with 
infrared point sources and image. They do not satisfy the 
luminosity criterion but have large line widths, consistent with 
their high masses. Both sources are located in the Cygnus X 
molecul ar cloud complex and were mapped in 1.2 mm con- 
tinuum (i Motte et ail 12007*). yielding masses 8M0 and 23 Mq, 
respectively (see their Table 1 and Fig. 13). Here we suggest 
that both cores are pre-UC H 11 regions. 

IRAS 21391+5802: '^CO emission detects a distinctive core 
and a belt of gas distributed along the southeast to northwest 
direction, which is coincident with the MIR background and 
several point sources. The core mass Mlte (1-3 x IO^Mq) is 
roughly half of its virial mass Mvh- (3.1 x IO^Mq), indicating 
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that the core is not yet gravitationally bound, responsible for a 
large line width (2.78 km s"'). Taking its relatively small size 
(0.32 pc) into account, we suggest that it is a star forming cluster 
where high-mass stars could eventually form. A 22 GHz water 
maser and outflow were identified in this region jPalla et al.L 
ll99lHZhang et al.Ll2005h . 

IRAS 22198+6336: MSX data are unavailable near this 
region, but the IRAS source matches the '^CO core peak well. 
A core as massive as 1.7 x IO-'Mq agrees with the conclusion 
deduced from the luminosity/line-width criteria. A 22 GHz wa- 
ter m a ser and outflow wer e identified in this region (Pal la et all 
119911: IZhang et all l2005h . We therefore suggest that it is a 
pre-UC H n region. 

IRAS 22506+5944: '^CO core coincides well with a lumi- 
nous IRAS point source and a bright MSX counterpart. A core 
as massive as 1.0 X 10^ Ms agrees with the conclusion deduced 
from the luminosity/line-width criteria. Thu s, we suggest that 
it is a pre-UC Hii region. Outflow activity (Wu et al , 2005b; 
Zhang etal., 2005) and a 22 GHz water maser (Palla et ay, 
199lh were identified, indicating that active star formation 
process is underway. 

In summary, '■'CO 7 = 2-1 mapping reveals at least 18 mas- 
sive cores from 14 maps. By means of individual analyses, we 
identify eight pre-UC H ii regions and one UC H ii region, two 
high-mass cores earlier than pre-UC H ii phase, four possible star 
forming clusters, and three sourceless cores. 
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Note. — Column (1) is source name shown in order of ascending IRAS number. Columns (2) and (3) list J2000 equatorial coordinates of each source. Column (4) lists the main beam temperature 
of '^CO y = 2- 1 if available, and Cols. (5)-(7) list the main beam temperature, local standard of rest velocity and line width of '^^CO 7 = 2-1, with Icr rms level in parenthesis. Column (8) is 
kinematic distance except for otherwise labeled. Column (9) lists the bolometric luminosity with asterisk (*) if one (or two) IRAS band flux is upper limit. Column (10) is excitation temperature. 
Column (11) is the opacity of '^^CO. Columns (12)-(13) pr esent average '^CO and H2 column densities. Column (14) l ists former works i f ava ilable. 

References. — (1) 6.7 GHz methanol maser survey dSzvmczak et all I2OOOII : (2) 22 GHz water maser survey dWouterloot et al.L [199311: (3) 22 GHz w ater maser jPalla et all Il99lh : (4) 
1665/67 MHz OH ma ser survey dWouterloot etail |1993|); (5) HCO+ 7= 1-0 survey (Richards et all 1 1987h : (6) CS 7 = 2-1 suryev dsronfman et all [1993): (7) NH^, (7,^) = (1,1) (2,2) survey 
I Molinarietal.Vl996'); (8) CO 7 = 2- 1 outflow mapping dZhang et aUl2005l : lKim & Kurtzi.i200d) : (9) CO. CS mapping jAo et alll2004h : (10) CO mapping fcuan et al.L [20081) : (11) CO mapping 
I Lee & Jung, 2003). N means non-detection. 

" f distance from " IWu et al.l ( l2004l) . '' IWouterloot et al.l ( Il993l) , '^ iMolinari et aLl ( ll996h , " iHariu et alj ( ll993l) . '= lMatthewsl ( [T97^ , iMotte et al.l ( l2007b : 

* kinematic distance unavailable, set as 1 kpc; 
mapped source; 

" more than one velocity component, only the strongest one listed; 
" upper limit. 



